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Abstract 

We  describe  the  speed  and  structure  of  steady- 
nonlinear  waves  moving  with  uniform  velocity  In  an  oblique 
direction  across  a  magnetic  field  in  a  collision-free 
plasma  with  temperature.   This  is  a  self -consistent  formu- 
lation in  which  the  plasma  motion  is  governed  by  the 
Vlasov  equations  for  ions  and  electrons  coupled  to  the 
Pols son -Maxwell  equations.   This  difficult  problem  can  be 
solved  in  the  limiting  case  of  small-amplitude  nonlinear 
waves  if  we  assume  that  the  same  scaling  laws  hold  for 
the  warm  plasma  as  those  for  a  cold  plasma. 
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1 .   Introduction 

The  effect  of  temperature  on  the  structure  of  nonlinear 
steady  plane  hydromagnetlc  waves  which  can  propagate  trans- 
verse to  the  magnetic  field  In  a  collision-free  plasma  has 
been  studied  In  [l].   This  Is  a  self -consistent  formulation 
In  which  the  equations  of  motion  (Vlasov  equations)  for  Ions 
and  electrons  are  coupled  with  the  Polsson  and  Maxwell 
equations.   In  the  limiting  case  of  a  wave  of  small  amplitude, 
a  particular  nonlinear  progresslng-wave  solution  In  the  form 
of  a  single  pulse  Is  obtained  by  using  the  same  scaling  laws 
which  hold  for  the  solitary  wave  In  a  cold  plasma  derived 
m  [2]. 

In  this  paper,  we  summarize  our  results  leading  to  a 
description  of  the  speed  and  structure  of  steady  nonlinear 
waves  moving  with  uniform  velocity  in  an  oblique  direction 
across  a  magnetic  field  in  a  collision-free  plasma  with 
finite  temperature.   Again  we  assume  that  the  same  scaling 
laws  hold  for  the  warm  plasma  as  those  for  a  cold  plasma; 
for  the  oblique  case  the  cold-plasma  scaling  laws  are  given 


1.  C.  S.  Gardner  and  G.  K.  Morlkawa,  "Effect  of  Tempera- 
ture on  the  Width  of  a  Small-amplitude  Solitary  Wave  in  a 
Collision-free  Plasma,"  Gomm.  Pure  and  Appl.  Math.  l8,  35-^9 
(1965).  — 

2.  G.  S.  Gardner  and  G.  K.  Morlkawa,  "Similarity  in 
the  Asymptotic  Behavior  of  Golllslon-free  Hydromagnetlc  Waves 
and  Water  Waves,"  New  York  University,  Courant  Institute  of 
Mathematical  Sciences,  Report  TID-6l8^  (May  i960). 


-  1  - 


in  [5].   Compared  to  the  transverse  case,  the  distribution 
functions  for  ions  and  electrons  have  a  considerably  greater 
variation  in  structure  in  the  oblique  case.   As  in  the  cold- 
plasma  oblique  case,  there  exists  a  critical  transition 
angle  at  which  the  steady  wave  shape  changes  from  a  positive 
(compression)  pulse  into  a  negative  (expansion)  pulse;  for 
the  warm  plasma  this  transition  angle  depends  on  the  tempera- 
tures of  both  the  ions  and  electrons  as  well  as  on  the  mass 
ratio  of  electron  to  Ion.   In  contrast  to  the  cold-plasma 
oblique  case,  singularities  in  the  higher-order  perturbed 
states  of  the  distribution  functions  arise  in  velocity  space. 
These  singularities  are  related  to  those  particles  which  move 
with  the  wave  in  the  direction  of  the  unperturbed  magnetic 
field.   Consequently  there  is  a  nonuniqueness  in  the  specifi- 
cation of  both  the  speed  and  structure  of  the  wave. 
Apparently  we  must  study  more  extensive  problems  which  include 
a  greater  variety  of  motions  in  order  to  completely  resolve 
this  difficulty,  e.g.,  time -dependent  problems  such  as  in 
[2]  and  [3]. 


5-   K.  W.  Morton,  "Finite  Amplitude  Compression  Waves 
in  a  Collision-free  Plasma",  Phys.  Fluids  7,  I8OO-I815 
(1964). 
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2 .   Scaling  Transformations  and  Equations  of  Motion 

We  consider  a  plane  plasma  wave  moving  with  uniform 
speed  Uq  in  the  negative  x  direction.   The  undisturbed 
magnetic  field  ahead  of  the  wave  has  constant  strength  Bp^ 
directed  at  an  angle  a  with  respect  to  the  wave  front;  a   is 
in  the  x-z  plane  (Pig.  l).   We  formulate  the  problem  in  the 
frame  of  reference  fixed  to  the  wave;  in  the  wave  frame 
there  is  a  constant  electric  field  in  the  y-direction  of 
magnitude  Eq  =  U^Bq.   Then  the  magnetic  and  electric  field 
components  may  be  written 

(B^,By,B^)  =  B^Csin  a,  b,  cos  a+B)  (l) 

and 

(E^'Ey^E^)  =  UqBq(E,  cos  a,  O)  (2) 

where  b,  B,  and  E  are  the  perturbed  field  components. 

We  look  for  nontrivial  solutions  which  go  to  a  constant 
state  at  X  =  -oo,  where  the  distribution  functions  f   for 
ions  and  f_  for  electrons  become  maxwellian  in  the  undis- 
turbed plasma  and  the  plasma  is  charge-neutral,  n_  -  n  =n^. 
The  overall  procedure  for  solving  this  problem  is  similar 
to  that  described  in  [ 1 ]  for  the  transverse-magnetic-field 
case.   In  the  oblique  case  we  also  make  a  coordinate  trans- 
formation in  velocity  space  to  facilitate  the  integration 
of  the  Vlasov  equations,  i.e.,  cylindrical  polar  coordinates 
(0,V,w)  in  which  the  w-axis  is  aligned  with  the  direction  of 
the  magnetic  field  B^  (see  Fig.  2): 
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u  =  U„(l+u  cos  a +w  sin  a)  (5-1) 

V  =  UqV  (3.2) 

w  =  U^(-u  sin  a  +  w  cos  a)  (3-3) 


where 


u  =  V  cos  0    ,  V  =  V  sin  9    . 

Guided  by  the  asymptotic  method  of  solution  In  [1], 
[2]  and  [3]  we  define  a  small  parameter  e  (related  to  the 
strength  of  the  wave)  In  terms  of  a  Mach  number: 

M^  =     g  Q  g  =  1  t  e^  (4.1) 

A  +  a 

where  A  Is  the  Alfven  speed, 

^2 

A^  =  5 {k.2) 

M.nQ(m_^  +  m_) 

and   a   Is  akin  to  the  sound  speed  and  must  be  determined 
as  part  of  the  solution.   In  the  transverse  magnetic  field 
limit,  3t  =  0,  we  must  get 

2  _  2(p^  +  p_) 
^  =  nQ(m^  +  mj  ^5) 

appropriate  for  a  plasma  with  specific  heat  ratio,  c  /c   =  2 

In  the  definition  of  e  In  (4.1)  we  consider  the  possibility 

2  - 
Of  having  both  subcrltlcal  (l  -e  )^  and  supercritical 

2  - 
(1  + e  )^  wave  speeds,  M. 
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We  now  assume  that  the  same  scaling  laws  and  expansion 
of  the  field  quantities  (with  respect  to  e)  hold  for  the 
warm  plasma  as  for  the  cold  plasma  [3]: 

where  (m  uVeBpi)  is  the  ion  gyro  radius  and  the  perturbed 
field  expansions  of  the  field  components  B  ,  B^ ,  and  E^, 
respectively,  are 

b  =  e\^^^   +   e^b^^)  +  O(e^)  (7-1) 

B=  B^B^^)  +  e^B^^)  +  O(e^)  (7-2) 

E  =  e^E^^)  +  e^E^^^  +  0(£^)  .  (7.3) 

The  superscript  quantities  are  independent  of  e .   We  intro- 
duce dimensionless  functions  f  and  f_  by 

f^  =  4  f^  and   f  =  -I  f  (8) 

We  need  only  work  with  the  distribution  function  f  for 
ions  since  the  Vlasov  equations  show  that  f  can  be  obtained 
from  f  by  replacing  terms  with  S(  )/'dx  in  f  by  (-m)S(  )/5x 
and  terms  with  /(  )dx  in  f_^  by  (-m~  )  /  (  )dx  where  m  =  m_/m_^ 
is  the  mass  ratio  of  electron  to  ion.  We  drop  the  subscript 
on  f  =  f  and  the  Vlasov  equation  of  motion  for  ions  can  be 
written 
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||-  =  e(l  +  w  sin  a  +  V  cos  9   cos  «)  ||  +  E^^  |^    (9-1) 


a       a^    (E.  -wb) 


-(e.  w^^(^bI--w)^ 


V 


j (f  cos  0) 


„rS  ,   .    /   ^    ^^  B  >,  ,  (l+w  sin  a)  ?/„   .   n\ 
^  r  (wb  -Ej 


+  Be 


V 


(f  sin  0)  -  Bf  cos  a 


^(^  +  Y^"  "^  (f  cos  0)! 


where  the  notation  E  and  Ej_  have  been  Introduced  to  simplify 
the  equation, 

E  =E   sin  a  +  b  cos  a   and   E,  =E  cos  a  -b  sin  a  .     (9-2) 

Recalling  the  expansion  (7)  we  see  that  the  right  hand  side 
of  (9.1)  is  0(e). 

With  the  velocity-space  transformation  (3)  and  scaling 
transformations  (3)  to  (6)  the  Poisson  equation  becomes 

00  27r  °° 


e(l  +  m) 


2       f    ^  ^ 

A  dE 

~7   dx 

c       J 

0  0 


(f^  -  f  )V  dV  d0  dw 


'  + 


(10.1) 


=  N(x) 


where  m  =  (m_/m  )  and  c  is  the  light  speed.   By  (l)  to  (6) 
the  Maxwell  equations  become 


6  - 


e(l+m)[l-(l+e^)  %]  g|  =  (l  +  £^) 

U 


+  w  cos  a)(f^  -  f_)V  dV  d0  dw 


00  2Tr  °° 
(      ( 


0   0    -00 


-V  COS  0  sin  a 


(10.2) 


=  (1  +  £^)j(x) 


and 


e(l  +  m)(l-(l+  e  )  %]  dx 

^0 


(l+e^) 


00    2^  °^ 

r     r     f 


J 

0   0    -00 


(V  sin  0)(f ,  -  f_)V  dV  d9  dw      (10.5) 


■  + 


J 


=  (1  +  e^)j(x)  . 

For  later  sections  we  use  the  abbreviated  notations: 

N  is  the  excess  number  density  of  singly  charged  ions  to 

electrons ; 
j   is  the  z -component  of  the  current  density;  and 
J   is  the  y-component  of  the  current  density. 
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3.    Series  Solution  of  the  Vlasov  (Collision-Free  Boltzmann) 
Equations 

We  assume  that  the  distribution  functions  f  =  f ,  and  f_ 
can  be  represented  as  power  series  expansions  in  e,    as  in 
[1]: 

f  =  fQ  +  8f^  +  E^fg  +  e^f^  +  £^^4  +  £^£"5  +  O(e^)     (11) 


and  similarly  for  f_.   f   is  prescribed  to  be  Maxwellian . 
If  now  we  had  followed  the  algorithm  for  determining  f,  in 
[1],  we  would  insert  (ll)  into  (9)  and  equate  coefficients  of 
like  powers  in  e;  then  we  would  l)  integrate  the  equation  for 
Sf,/S9  beginning  with  k  =  1,  and  2)  use  a  compatibility  condi- 
tion for  the  solvability  of  the  next-order  equation  for  f,^, -,  , 
i.e.,  integrate  the  equation  for  Sf,  -,/S0  from  0  =  0  to 
0   =   27r,  in  order  to  evaluate  the  arbitrary  ©-integration 
variable  for  f,  in  l).   This  procedure  is  correct  in  principle 
but  unnecessarily  tedious.   A  much  simpler,  but  equivalent, 
algorithm  for  determining  f,  is  available  by  recognizing  the 
unique  functional  form  of  f  with  respect  to  0,    namely, 

k-1 
^k  =  \o  +  ^  ^\£   ^°^  ^^  +  \£    ^^"^   ^^^  ^12) 

i/^—  -J- 


and  obtaining  (from  the  Vlasov  equation)  recurrence  relations 

relating  the  coefficients  A,  „,   B,  „,    and  A,  ^,  which  are  func- 

ki'      k£  kO 

tions  of  V,  w,  the  perturbed  fields  and  their  derivatives,  to 
the  lower  order  coefficients. 
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We  prescribe  the  unperturbed  distribution  function  f^ 
to  be  maxwellian 


^0  =-   ^00  =  5(V^  +  w^)  (13.1) 

where 

m  U^ 
?  =  $+  =  (yA)^/2  exp[-Y(v2+w2)]  ,  y^   .  ^=  Y  (13.2) 


For  electrons  the  unperturbed  distribution  function  is 

$'_  =  (t./tt)^/^  exp[-7_(v2  +  w2)]  ,   T_  =  ^  •     (13.3) 


We  note  that  a  more  general  anisotropic  distribution  f  , 

in  place  of  (13-1),  is  allowed,  i.e.,  a  monotonically 

2     2 
decreasing  function  of  (v  +  Aw  ),  A  >  0.   The  perturbed  f  j^ '  s 

(k  =  1,2,3,4)  which  we  need  in  the  later  sections  are  the 

following: 


fg  =  A^Q  +  Ag^cos  e  (15.1) 

X 

A^Q--  ^  (B^^Vcos  a+2  ^(^^w),  I^' "^  =  I  E^^-^'^dx   (15.2) 


-00 


Ag^  =   2B^^)  $''V  (15.3) 

where  W  =  (l  +  w  sin  a)  and  primes  indicate  derivatives  and 
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f      =   B^^sln    9   +  3^2^^"  ^^'      ^^£   =    °'      -^   =   0,1,2,...  (l6.l) 


2        2 

(l)>T.',w^,,       V  COS  a 


B^n    =  b'  ^•^''^V(2W 


^31 


W 


)+ -r (16.2) 


W 


B,p   =^B'(^)^'v2 


^2         2 


cos  a 


(16.3) 


£"4   =   -^40  "^   Ai^^cos    9   +   A|^2^°^   ^0   +    A^^^cos    39 


(17.1) 


Bi^^   =    0    ,       ^    =    0,1,2,3,.  . 


T^"  (1)   (Ty'T.2              /,,        3  V^cos^a 
A40  =  B  I  V   cos   a(W    -  ^  ^ 

(1)     _         '(1)      s     .T,'..2 


(Ej_^^    -   b    ^^w)   ^  V   cos  a 


+ 


W 


+    (B(l))'{fV2(l   + 


(17.2) 


2         2 
V   cos   g 

2 
2W 


.TVrw        ^'^^    -I  Sing)        ^^^^^2      ^^^2 
+  $*[¥  +   ^^ +  ^ ] 


W 


4W 


3 


+ 


,    (1)    M)   V^cos   g   j_    .r. 


X 


=?'f<^W)ii,(£: 


_  W\ 


+   2( 
-  i^    (B^^  Vcos  g   +   2   l^^K) 
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A^^  =  -b"(1)  f'v(2W^-  Iv^oos^a)  -2(E;tl)-  b'(l'w)^V 

-  (B(l))2vcosa(4$V24j^.i^^%^)         (17.3) 

w 


-I  (e;(1)-  b'(l)w)  i:x!£2^  (17. i^) 


^43  ^  "  T2  ^"^^^   ^'V^cos^a  (17-5) 


By  Inspection,  (13)  to  (l?)  show  that  f^  can  be  distinguished 
between  even  and  odd  k.   For  even  k 

2k-l 
f2k  =  hl.,0   +  ^  ^2k,^  ^°^  ^^  '   ^0  =  ^00        (lS-1) 

and  for  odd  k 

2k 


f2k+l  =  ^'^2k+l,^  ^^"^^  ■        fl  =  °  (18.2) 


These  particular  forms  (l8)  for  f,  are  useful  In  the  next 
section  in  determining  the  expansions  of  the  moments  of  f 
and  f   in  powers  of  e .   We  point  out  that  the  singular 
integrals  which  arise  when  we  calculate  the  moments  of  f 
and  f  come  from  the  zeros  of  powers  of  W  =  (l  +  w  sin  a) 
occurring  in  the  denominator. 
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4 .    Computation  of  the  Moments:  Number  Density  N  and 
y-  and  z -Components  of  Current  Density  J  and  j 

Because  of  the  forms  (l8)  of  the  perturbed  distribution 
function  f,  for  even  and  odd  k,  the  moments  N,  j,  and  J  In 
(10)  have  the  following  expansions  with  respect  to  e: 
Only  fp,  contribute  to  the  number  density  N, 


N  =  e^N^^)  +  e\^^^  +  O(e^)  (l9.l) 


where  N^  '^  =  0  since  the  undisturbed  plasma  (x  -    -co)  Is 
charge -neutral ;  and  also  only  fp,  contribute  to  the  z- 
component  of  current  density  j, 


J  =  .2J(l)  +  s^•(2)^-o(.6).  (19.2) 

Only  fpk+l  "^o^trlbute  to  the  y-component  of  current  density 
J, 

J  .  e^j(^^  +  £5j(2)  +  o(e'^)  .  (19-3) 

By  the  perturbed  field  expansions  (7)  and  the  perturbed 
density  expansions  (19)^  the  Pols son -Maxwell  equations  (lO) 
become 

s\i  +  „,)  4^=  e2(N(lUA(2))+o(E«)       (20.1) 
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e\l+m)[l-(l+e2)  ^]  ^^  (20.2) 

^0 


,  e^(l  +  ,2)(j(l)  +  e2j(2))  +  0(efi) 


2 
I 

U 


2 
e^(l  +  m)[l-(l+e2)  ^]  ^  (b^^^  +  e^B^^))         (20.3) 


0 
=  e^(l  +  e^)(j^^)  +  e^J^^h  +  O(e^) 


Now  the  perturbed  quantities  with  superscripts  are  independent 
of  e  and  we  can  equate  coefficients  of  equal  powers  in  e  in 
(20).   To  lowest  order  in  e  we  get  certain  consistency  relations 
including  the  speed  of  a  weak  steady  wave  in  the  limit  e  — >  0. 
To  the  next  order  in  e  we  get  differential  equations  describing 
the  wave  structure. 
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5.   Consistency  Relations  and  Speed  of  a  Weak  Wave 

To  the  lowest  order  In  e  the  Pols son -Maxwell  equations 


(20)  yield  to  0(e  ), 
n(i)  ^  0 


(21.1) 


j(^)=o 


(22.1) 


and  to  0(e), 


a 


(l  +  m)(l  -  \)  B 
Uo 


(1)  .  j(l) 


(23.1) 


By  (10.1)  and  (I9.I),  the  first  order  number  density 


N^^)  is 


N' 


■^)  =   ff  (f^  -  f2)V  dV  d0  dw 
=  \\l    (A+o  -  A^q)  V  dV  de  dw 


(21.2) 


Carrying  out  the  Integration  with  A^^^  given  in  (15),  (21) 
states 


20 


B^^-'sin  a  cos  a(j6't  -j6^  )  +  2^*'^^[y^(1  -jt)'t' 

1        1  -^11    +     1 


(21.3) 


-1 


+  m  ^Y_(l-/5^)]  =  0 

where  ^-.    is  the  principal  part  of  a  singular  Integral 

2 


-CD 


e 

(1  +  KO)) 


-1  .  2 


^1  -  r'^P]      TT^i;^  do)  ,    k'"  =  Y"^sin^a  (24.1) 


K    ^^P 


(-l/K^)  I 


■l/K 


exp(a  )da 


0 
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For  p   order  singularities  we  define  fi     as 


-co 


(1  +Ka))P 


doo 


(24.2) 


is 


By  (10.2)  and  (19-2)  the  first  order  current  density  j 


(1) 


(1) 


=  jj I  (-V  cos  0  sin  a  + w  cos  a)(f2  -  fg)^  dV  d9  dw 
^  \\\    ^"2"  4l  ^  ^^^   ct  +A+Q  w  cos  a)  (22.2) 

-  (-  ^  k~     V  sin  a+  A2Q  w  cos  a)  ]V  dV  dS  dw  . 


Carrying  out  the  integration,  (22)  states 


(1) 


+ 


B  '^  '^sin  a  cos  a  (/)^  -  jb^) 
.(1), 


(22.5) 


+  2eJ^''[y_^(1  -  ^\)   +   m-^Jl  -  ^l)] 


=  0 


which  is  identical  to  the  x-derivative  of  (21.3)-   Thus 

j^-'-''  =  0  is  consistent  with  N^  -^  =  0  but  gives  no  new  information 

(1) 


By  (10.3)  and  (19.3)  the  first  order  current  density  J 


IS 


j(l)  =  11   (V  sin  e){f't    -   f^)V  dV  d0  dw 


(23.2) 


\\\  <b; 


^.V"^  -  B,.V^)dV  d0  dw 
31      ^1 


Carrying  out  the  integration  with  B^-j^  given  in  (16),  J^  ''  i£ 
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^ 


+ 


p'-. 


j(l)  =  -b'  ^^^[  (1  +m)  -  (-i  +  m  -^)  cos^a] 


T 


+ 


Y_. 


(23.3) 


IT      -L      -L 


sm  a 

Since  (21.3)  gives  a  linear  relation  between  E^  '  and  B  ^  ^ 
(23.1)  and  (23.3)  yield  a  relationship  independent  of  per- 
turbed field  quantities: 


2 


(1  +  m)  ^=  i(^+  m  J^) 


U 


0 


"T 


+ 


y-' 


(25) 


^+ 


-I 


O^l  -  J^IY 


Ly^(1  -;?5+)  +m"\_(l  -/^-)J 


2 
cos  a 


This  consistency  condition  determines  the  "sound"  speed   a 

in  (4)  and  relates  it  to  the  other  physical  parameters  a,  m, 

U„,  T_,  and  T  .   In  the  limit  e — >0  in  (4),  the  weak  wave 

speed  is  given  by 
,2 


^0 


(26) 


combined  with  (25).   Eqs .  (25)  and  (26)  recover  the  results 
in  [1],  [2],  and  [5]  for  the  limiting  cases  of  l)  transverse 


magnetic  field,  a  =  0,  and  2)  the  cold  plasma,  T_ 


+ 


0, 
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6.    Differential  Equations  for  the  First  Order  Field 
(e(^),  b(^),  B^^)) 

From  the  Polsson-Maxwell  equations  (20)  the  system  of 
equations  which  we  get  to  the  next  significant  order  In  e 
are  three  differential  equations  for  six  unknown  field 
quantities  (E^-^S  E^^K    b^^S  b^^^,  B^^S  B^^^).   But  because 

of  the  asymptotic  nature  of  our  approximation  procedure  the 

(2)         (2)         (2) 
second  order  fields  (E^  S  b^  ',    B^  ^)    can  be  eliminated  by 

the  consistency  relations  (21.5)  and  (25).   Hence  we  obtain 

three  differential  equations  for  only  the  first  order  fields 

(E^  S  b^  S  B^  ^).      This  system  of  three  equations  can  be 

reduced  further  to  a  single  differential  equation  for  one 

field  component,  say  B^   .   Since  the  complete  equations  are 

quite  lengthy  we  will  only  show  the  essential  elements  In  the 

following  description  and  not  the  explicit  details. 

The  Polsson  equation  (20.1)  gives  to  0(e  ): 

(1+m)  4e'^^^  =  N^2)  (27.1) 


2 
c 


where 


N 


^^^  =    I  (^t  "  ^4^^  ^^  ^^  ^^  (27-2) 

J  J  J 

=    \j\     (4o  "  ^40^^  ^^  '^^    ^"^ 

=  q,B"^l)  +^  (B^l))2  4-  q^E'^^)  +  q^b'^^) 

+  m-\    (1-/5:)]  ^£5^]  , 
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whe 


re  q.,  ,  qp ,  q^,  and  q^  have  known  dependence  on  the  pre- 


scribed physical  parameters  (a,  m,  U„,  T_,  T  );  In  particular, 


q^  =  0, 


The  Maxwell  equation  (20.2)  for  the  z-component  of  the 


current  density  j  gives  to  0(£  ): 


(28.1) 


where 


.:(2) 


(-V  cos  e    sin  a  +w  sin  a)  (fj  -  f  ^^  )V  dV  d0  dw    (28.2) 

fff     ^41  + 

[  (-  —2—  V  sin  a  +  Kj^q   w  cos  a) 

-  ( ^  V  sin  a  +  ku^   w  sin  a)]V  dV  d0  dw 


=  r^B 


2 


sm 


+  m-^Y  (1  -Jd-)]   22^^] 

where  r-,  ,  r,, ,  r^,  and  r^  depend  on  the  given  physical 
parameters . 

The  Maxwell  equation  (20.3)  for  the  y-component  of  the 
current  density  J  gives  to  0(s  ): 


2 
(1  +  ni)[l  -  ^)  b'^^)  7  B'^^)] 

"o 


J 


:2) 


(29.1) 


where  from  (4)  the  (-)  sign  refers  to  supercritical  wave  speed 
and  the  (+ )  sign  refers  to  subcrltlcal  wave  speed;  and 


-  1^ 


j(2)  =  jjl  (V  sin  0)(f^  -  f^)V  dV  d0  dw  (29-2) 


J  J  (V  cos  0)(^  -  ^)  V  dV  d0  dw 

-  [A^^VW  +  (A^Q  +  -2^)  V^cos  a]  I  V  dV  d0  dw 

-  B^^^J^"^^  cos  a 


S3_b'"(1)  +  S2b(1)b'(1)  +  s/(l)  +  s^b"(l) 


+  {B'^2)[-(l+m)  +  (-1  +  m  J^)  cos^a] 

-  e(2)(^+  ^^-)  cos_^  ) 
II   ^'^1   ^1^  sin  a,  } 

where  s-,,  Sp ,  s^,  and  s^  depend  on  the  given  physical 
parameters.   In  the  calculation  of  J^  '    the  explicit  need 
for  fp-  was  circumvented  by  taking  the  (V  cos  9) -moment  of 
the  Vlasov  equation  (9-1)  and  using  the  following  moment 
relation  (obtained  by  an  integration  by  parts): 

I J j  (V  cos  0)  |§  V  dV  d0  dw 

=1  (V  sin  9)fV  dV  d0  dw  .      (29-3) 

As  shown  by  (29.2),  {'^f  /'bQ)    can  be  expressed  in  terms  of  f, 
with  k  <  4. 

The  elimination  of  second  order  quantities  (e^  S  b '  ' , 
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B^^^)  from  (27),  (28),  and  (29)  is  easily  described, 

recalling  the  definitions  of  E  and  $  from  (9-2)  and 

It       11 

(15.2):   1)   Clearly,  the  second  order  terms  can  be  elim- 
inated between  (27)  and  (28)  since  the  expressions  in  curly 

brackets  are  identical;  2)  taking  the  derivative  of  (27) 

(2  ) 
and  eliminating  E^  '  using  (29),  the  coefficient  of  the 

(2) 
remaining  second  order  term  in  B '     becomes  zero  by  the 

consistency  condition  (25);   and  3)  the  third  Independent 

equation  for  (E^  S  b  ^  '' ,  B^^)  is  (22.3).   From  these  three 

differential  equations  we  can  obtain  a  single  third-order 

nonlinear  differential  equation  satisfied  by  B^  ^    of  the 

form 

Rb"'^^)(x)  +  (;  1  +  SB^^^x))b'^^)(x)  =  0         (30.1) 

where  the  coefficients  R  and  S  are  given  by 

(1  +m)R  =  (30.2) 

=    (l+m5)+^    (^+^)(_ll+23    sln2a)+|   (4+   m^  4)^°^^^ 

1  2    - 

p  p^-'^i'^^i  ?T~'+         -        2 

+    Pcos  a[-(l  -m^)  +  ^   (^  -  -y-i)    cos  a +\-{j6^  +  m/)^)oos  a] 

-[(1  -m2  +  l(^   -  ^)(-3+5   sln^a)  -  i(l  +  m)  Tcos^aJ^ 

+ ^ 1 

cos^aUl  +  m)  +  [(4  +  ^)   +  ^   K  -^l)]sln^al 


20    - 


(l+m)S  =  (50.5) 


^+ 


cos  ab(l  +  m)-[(;^  +  -;^)    +  V   ('^1  ->^i)]c°s  a 

-  (t^  +  -^)  (2-5   sin  a)   +   3(— -  +  m  — )    cos  a 

T+  r_  1-+  r_ 

_  I    (^  +  !1^)    sm^a    cos^a  -  V  {P\  -  P\)  {-^  -^  sin^a) 
+         -  ?  "P      -^^^      '^^^  2  2 

+  r(i^2  -^^2)''°^  ^  +  2  ^Y  "y  )  ^^"  ^  "^^^  ^ 

-  ^  [(i2^2   +  "^"^^2)  -5(^3  +  m"^j25^)]sln2a   cos^a 

-  \-  [(7/2    "   =^~^Y_/^2)  -  (7/3  -m~2Y_/~)sln2a    cos^aj 


Where  p      Is  defined  by  (24.2)  and  specifically  p,^   =   -2k   (l  -p^) 
and  /!)^  =  -Yr^{p-^  -^^)    and 

P=  '^   ,^  (30.4) 

[y^(1  -;25+)+m-lY_(l-j2)-)] 


A  preliminary  study  based  on  suggestions  by  C .  S. 
Gardner  (private  communication)  indicate  that  using  these 
relations  for  j6^   and  jt)^   together  with  the  principal  value 
integral  j25.  in  (24.1)  is  equivalent  to  taking  "cut-off"  dis- 
tribution functions  (at  x  =  ±00)  which  eliminate  particles 
which  turn  around  or  are  trapped  by  the  steady  wave. 
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The  differential  equation  (30)  Is  of  the  same  form  as  that 

derived  In  [l]  and  [2]  for   l)  the  steady  wave  moving  In  a 

transverse  magnetic  field  and   2)  the  cold-plasma  oblique 

case  [3].   The  (jjl)  signs  allow  for  both  super-  and  sub- 

crltlcal  wave  speeds.   Although  the  coefficients  R  and  S  are 

lengthy,  the  transverse  limit  a  =  0  and  the  cold-plasma  limit 

T  =  T  =0  are  easily  checked  by  inspection. 

We  have  obtained  the  parametric  dependence  of  R,  S,    and 

the  wave  speed  (25)  on  a ,  T_,  and  T  by  electronic  computation 

for  three  interesting  cases:   T_  =  mT   (cool  electrons), 

T  =  T,  (equal  temperature),  and  T^  =  0  (cold  ions)  (see  Pigs. 

3,  ^,    5,    and  6).   For  T_  =  mT  ,  we  have  Y_  =  Y,  =  J, 

j6~  =  j6'^  =  j6     and  P=  0;  for  T_  ^   T  ,  we  have  y_  =  my^^ .   In 

both  of  these  cases  the  critical  transition  angle  a   where 

°    cr. 

the  solution  changes  from  a  positive  pulse  moving  with  super- 
critical velocity  to  a  negative  pulse  traveling  at  subcrltical 
speed  Increases  rapidly  as  a  function  of  temperature.   However, 
the  cold  ion  case  (T   =  0,  T  >  O)  yields  a  singularly  different 
result:  a    is  practically  independent  of  T  ,  the  transition 
angle  being  approximately  the  cold  plasma  value  found  by 
Morton  [3]-   The  strong  dependence  of  the  coefficient  R  on 
the  ion  temperature  T,  is  emphasized  in  Pig.  4;  the  width  of 
the  pulse  is  proportional  to  R^ .   The  coefficient  S,  which  is 
proportional  to  the  amplitude  of  the  pulse,  exhibits  some  non- 
uniform behavior  for  both  small  and  large  obliqueness  angle  a. 
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particularly  with  appreciable  electron  temperature  T   (Fig.  5) 
This  non-uniform  behavior  of  S  Is  also  reflected  In  the  wave 
speed  (U  /a)  (Fig.  6).   However,  for  the  cool  electron  case 
(T_  =  mT  )  the  wave  speed  transverse  to  the  undisturbed 
magnetic  field  (U  cos  a/a)  Is  practically  Independent  of  the 
obliqueness  angle  ct  over  the  calculated  range  of  Ion  tempera- 
tures . 


7 •   Pressure  Tensor ,  P 


ij 


The 


pressure  tensor  for  Ions  (or  P.  .  for  electrons)  Is 


where 


00 


J 


f  du  dv  dw 


-00 

00  27r  "^ 


!  1 1  -' 


f  V  dv  d0  dw 


O  O   -00 


u    uv  uw 

2 

q. .  =    vu   V  vw 

wu   wv  w 


\ 


/ 


and  In  cylindrical  polar  coordinates. 


(31.1) 


(31.2) 


u  =  V  cos  0    ,        V  =  V  sin 


(31.3) 
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u   = 


2  2 

^  (1  +  cos  29)  ,   v^  =  ^  (1  -cos  29) 


(31.4) 


uv 


V 


sin  29,   uw  =  wV  cos  9,   vw  =  wV  sin  9     (31-5) 


is 


The  power  series  expansion  In  e  for  the  pressure  tensor 


P.  .  =  P-&.  .  +  e^pf^)  +  e\{V    +   e^?{V   +   0(£5)     (32.1) 
ij      ij      ij       ij       ij      ^   ^     ^     ^ 


where  P  Is  the  dlmenslonless  scalar  pressure 


P  - 


(V^  +  w^ )$  V  dV  d9  dw 


(32.2) 


=  kT/m^ug 


1/2T 


+ 


(cf.  eq.  (13-2))  and  the  perturbed  pressure  tensor  Is  calcu- 
lated using  the  series  solution  for  f  given  by  f.  In  (13)  to 
(l8)  and  calculating  the  moments  (31)-   We  summarize  the 
results  for  P:.'^:  P:.'^  Is  given  by  P:.'^  =0  for  1  /  j  and  the 
diagonal  elements  are 


(2)  _  p(2) 


11 


=  P 


22     2 


HI 


B^^Vcos  a   ^,^^^(1  -jt).  ) 
± +    _. — 


A^QV^dV  d9  dw 


(1) 


(33.1) 


Y 


sm  a 
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,(2) 
33 


1 


A2qW  V  dV  d0  dw 


(33.2) 


^^  [B^"^^(/)-j^  -l)sin  a  c 


OS  a  + 


^^  ^sin  a 


sm'a 


+  2y^|_i)(i  -jb^)] 


X2)  , 


(1) 


(recalling 


P.^^'  is  directly  proportional  to  the  field  B 

(21.3))  and  comparable  to  the  usual  guiding  center  theory. 

v[^)    has  both  right  and  left  diagonal  elements  equal  to  zero. 


,  ^y.'    =  0  for  i  =  j  and/ or  (l  +  j)  =  ^  and  the  non- 


1  .e 

diagonal  elements  are 


>(3) 
12 


(3)  _  1 


21 


^ 


II 


B^oV  dV  dO  dw 
32 


(3^.1) 


1  B' 


(i: 


cos  a 


T 


,(3)  _  p(3) 


23 


32 


MK 


B^.wV  dV  d0  dw 
31 


(34.2) 


B 


2  Y  sin  a 


1)0  2 

[sin  a  +  (yT-  2)cos  a] 


V).  is  directly  proportional  to  the  field  gradient  dB^  ydx. 
p'.''  has  non-zero  right  and  left  diagonal  elements  and  V)  .'  = 
for  i  /  j  and  (i+j)  /  4;  P:.^  contains  second  order  field 

;2) 


0 


quantities  (E 


:2)   k(2) 


B 


and  is  rather  lengthy  so  only 


the  integral  form  of  the  components  are  given 


-  25  - 


^11     ^22     ^ 


(2A^qV^  +  A^2^^)dV  de  dw       (35.1) 


?^V   =   ||[  A^^qW^V  dV  d0  dw  (35.2) 

^iV    ^  ^iV    ^  \  \\[   Ai^iWV^dVdSdw  (35.3) 


4 
Thus  to  0(e  )  the  pressure  tensor  P. .  is  filled  with 

symmetrical  non-zero  elements. 


8.   Comments  on  Singularities  In  Velocity  Space 

We  can  discuss  the  singularities  In  velocity  space  only 
briefly  here.   We  begin  the  discussion  by  assuming  that  the 
equilibrium  distribution  function  f^   and  the  first-order  dis- 
tribution function  f-,  are  spatially  uniform  so  that  there  are 
no  charge  nor  current  densities  In  the  equilibrium  plasma. 
Then, 

fo  -   Aqq  =  $(V^  +  w^)  (36) 

fl  =  A^Q  =  0  (37) 

f2  =  Agrees  0  +  A^Q  (38) 

=  2B^^^  ^  '  V  cos  0  +  Xg^^''^'^) 
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where  Ap„  =  J^p  is  the  9-lntegratlon  variable.   In  order  to 
determine  ^„   we  impose  the  condition  that  to  next  order  f^ 


is  continuous  in  9,    i.e.,  using  the  Vlasov  equation  (9)  and 

27r 

^  d9  =  0  (39) 


r  '' 


0 
and  we  get 


I 


(1  +  w  sin  a)   ^^^  =  -  B  '  ^  "^  ^$'v^cos  a  -2E^'^^$'w  .    ('10. l) 


Integrating  'Xp  using  the  boundary  condition  fp  =  0  at 
X  -    -00  we  get 


(1+w  sin  a)X2  =  -B^  ^ -^  $' 'v^cos  a  -  2|^^)$'w  .    (40.2) 

Since  the  range  of  integration  for  w,  when  taking  moments, 
is  -00  <  w  <  00,  the  quantity  W  =  (l  +  w  sin  a)  can  be  zero. 
Hence  we  Introduce  a  5-function  in  velocity  space  and  write 


Xg  as 


ApQ  =  X2  =  -  w  (B^^V^cos  a  +  2^^^K) 


+   B^-^^(x)-C(V,w)  •5(W)  (41) 


for  a  >  0  and  C(V,w)  is  an  undetermined  regular  function. 
The  difference  between  A^^   in  (4l)  and  ApQ  in  (15.2)  is  the 
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6-function  term.   We  take  the  6-function  term  to  be  propor- 
tional to  the  perturbed  field  quantity  B  ^  ;    otherwise  the 
speed  of  the  wave  {\]^  a),    cf.  (25),  will  be  x-dependent.   But 
now  the  algorithm  for  the  asymptotic  solution  given  in  the 
previous  sections  does  not  determine  C(V,w)  nor  two  additional 

regular  functions  of  (V,w)  arising  from  the  higher  order 

2  3 

singularities  when  W  =  0  and  W  =  0  in  f ^ ,  cf.  (l?)-   Thus 

there  is  a  parametric  non-uniqueness  in  the  determination  of 

both  the  speed  and  structure  of  the  steady  wave  in  terms  of 

the  given  physical  parameters  (a,  m,  T_,  T  ),  i.e.,  in  the 

determination  of  (Up/a)  and  the  coefficients  R  and  S  of  the 

third  order  differential  equation  (30).   Although  the  present 

calculations  apparently  can  be  interpreted  in  terms  of  cut-off 

(or  cut-out)  distribution  functions,  a  more  precise  analysis 

seems  desirable  in  the  future--and  possible  (see  footnote, 

page  21)  . 
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Figure  2o.  Velocity  Space: 
Rotation  by  Qt 


w  out  of  paper 


u 


Figure  2b.    VelOCity    SpOCe. 

Transverse  Plane  to  w 
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FIGURE  3.    CRITICAL    OBLIQUENESS   ANGLE  vs.  PRESSURE 
FOR    COEFFICIENT     R-0. 
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FIGURE  4.  COEFFICIENT  R   vs.  OBLIQUENESS   ANGLE   a 
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FIGURE  5.  COEFFICIENT   S    vs.  OBLIQUENESS  ANGLE   a 
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FIGURE  6.  (WAVE  SPEED)"    vs.  OBLIQUENESS  ANGLE  a 
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This  report  was  prepared  as  an  accoiuit  of 

Government  sponsored  work.   Neither  the 
United  State?,  nor  the  Commission,  nor  any 
person  acting  on  behall"  of  the  Commission: 

A.  Makes  any  warranty  or  representation, 
express  or  implied,  with  respect  to  the 
accuracy,  completeness,  or  usefulness  of 
the  information  contained  in  this  report, 
or  that  the  use  of  any  information, 
apparatus,  method,  or  process  disclosed 
in  this  report  may  not  infringe  privately 
owned  rights;  or 

B.  Assumes  any  liabilities  with  respect  to 
the  use  of,  or  for  damaces  resultinj';  from 
the  use  of  any  information,  apparatus, 
method,  or  process  dlf>clor>ed  in  ihU: 
report. 

As  used  In  the  above,  "person  acting  on  behalf 
of  the  Commission"  Include.",  any  employee  or 
contractor  of  the  Commission,  or  employee  of 
such  contractor,  to  the  extent  that  such  em- 
ployee or  contractor  of  the  Commission,  or 
employee  of  such  contractor  prepares,  dis- 
seminates, or  provides  access  to,  any  infor- 
mation pursuant  to  his  employment  or  c.ontrafrl. 
with  the  Commission,  or  hl;>  employment  with 
such  contractor. 
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